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In order to quantify entanglement between two parts of a quantum system, one of 
the most used estimator is the Von Neumann entropy. Unfortunately, computing this 
quantity for large interacting quantum spin systems remains an open issue. Faced with 
g this difficulty, other estimators have been proposed to measure entanglement efficiently, 

Q mostly by using simulations in the valence-bond basis. We review the different proposals 

^ and try to clarify the connections between their geometric definitions and proper ob- 

I I servables. We illustrate this analysis with new results of entanglement properties of spin 

1 chains. 

^ Keywords: entanglement measures; quantum spin models 
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•/^ 1. Introduction 

• Entanglement is a fundamental notion of quantum mechanics, that has over the 

recent years gained popularity as a way to provide new insights in the quantum 
many-body problem. From the condensed matter point of view, one of the most 
" ^ I interesting promises of the study of entanglement properties is the possibility to 

^ automatically detect the nature of quantum phases and of quantum phase transi- 

^> tions. In this approach, there is no need to provide a priori physical information 

or input, such as the specification of an order parameter. The detection can oc- 
cur through the study of the scaling (with system size) of various entanglement 
estimators. For instance, the scaling of the von Neumann entanglement entropy 
for one-dimensional systems is different for critical and gapped systems - allowing 
their distinction. For a recent review of various properties of entanglement entropy 
in condenser matter, see Ref. T, 

Unfortunately, computing von Neumann entropy for large interacting systems 
remains a highly difficult task, except in one dimension where it can easily be ac- 
cessed with Density Matrix Renormalization Group (DMRG) algorithm. As a conse- 
quence, several other estimators have been proposed recently, that can be computed 
efhciently with large-scale quantum Monte Carlo (QMC) techniques. Most of these 
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proposals use simulations in the valence-bond (VB) basis. In section [2j we define 
useful properties of the VB basis. In section |3j we review the different estimators 
for entanglement entropy and clarify their definitions, thus allowing to compute 
some of them with other techniques. Finally, in section [4] we present results for 
frustrated spin- 1/2 and spin-1 chains. 

2. Valence-Bond Basis 

Let us first consider the antiferromagnetic (AF) spin-1/2 Heisenberg Hamiltonian 

H = ^ Jij Si ■ , 

ij 

which conserves the total spin St of the system. AF interactions > favor a 
singlet St = ground-state. 

It is well-known that any singlet state can be expressed in the VB basis, where 
spins couple pairwise in singlets (|t|) — ||t))/^/2- The VB basis is overcomplete 
(there are more VB coverings than the total number of singlets). Another basis is 
the bipartite VB basi^, where the system is decomposed into two sets such that 
two spins forming a singlet necessarily belong to different sets. This basis is smaller, 
but still overcomplete. Choosing a bipartition of the A^-sites lattice into two equal 
sized subsets A and B (for instance different sublattices for a bipartite lattice) , then 
the bipartite VB subspace is generated by all the bipartite VB states 

\vv)^ (g) (2) 

where [i,j] is a SU(2) dimer state and I? is a dimer covering of the system. 

A crucial property of the VB basis is its non-orthogonality. Indeed, for any pair 
of VB states \(px>) and \(pv}, the overlap is given by 

(^P,|^p)^e(I?,I?')2"'(^•^')-^/^ (3) 

where n; (T>, V) is the number of loops of the overlap diagram obtained by super- 
imposing the two dimer coverings V and V and s {V, T)') a sign directly related 
to the antisymmetric structure of the two-site singlet. Since all sites of the overlap 
diagram are involved in exactly one dimer from T) and one dimer from T)' only 
even sized closed loops are possible. Each loop produces a multiplicative factor 2 
to Eq. ([3]) since only the two sequences f, i, • • • , i and |, t . . . , t along the loop lead 
to a non-vanishing contribution. For arbitrary range VB states, the sign 8(1), V) 
cannot be fixed to a constant independent of the couple [V, V) by a generic conven- 
tion. However, considering only bipartite VB states, a convenient choice consists in 
orienting the dimers from sublattice A to B. In this case, it can be readilly checked 
that e (I?,I?') = +1 hence making {ip-D'ly^v) an unsigned quantity which turns out 
to be crucial for Monte-Carlo applications. 
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Despite being an overcomplete and non-orthogonal basis, Sandvik has shown 
how to make large-scale QMC simulations in the VB basis'^ , thus allowing to sample 
ground-state of non-frustrated spin models. The basic idea consists in applying the 
hamiltonian a large number of times on a trial state, in order to project out onto the 
ground-state of the system. Moreover, recently Sandvik and Evertz have improved 
this algorithm by implementing non-local moves.'^ 

3. Measuring Entanglement 

3.1. Generalities 

Over the past few years, various measures of entanglement have been used to in- 
vestigate quantum phase transitions and states of matteiEEEIll The most famous 
estimator is the so-called von Neumann entanglement entropy (5''^^, see below for 
its definition) and in our context, one wants to compute how this entanglement 
entropy (EE) scales with the block size. While many interesting properties of EE 
have been derived exactly for integrable models or systems with exactly known 
ground-states (GS), the calculation of EE for general interacting quantum systems 
is an exacting task. On the practical side, even numerical simulations are difficult. 
Exact diagonalization (ED) is limited to small systems and cannot precisely verify 
scaling properties and EE is not accessible to QMC methods. In Id, EE can be cal- 
culated within the DMRG mcthocF' and crucial information on the criticality can 
be extracted using conformal field theory tools. ^ Unfortunately, DMRG is not avail- 
able in larger dimension. Therefore, several other estimators have been proposed to 
measure entanglement using QMC techniques in the VB basis. 

From the experimental point of view, if one considers small blocks (1 or 2 sites 
for instance), then the reduced density matrix (and thus EE) can be expressed in 
terms of all correlations on the same block: local magnetization, nearest-neighbor 
correlations ... so that in principle, local measurements can give access to these 
entropies. But, since we are mostly interested in large blocks to get insights on 
the nature of the quantum state, it seems impossible to be able to measure EE 
in this limit. Despite this negative result, Klich and Levitov have shown that, for 
free fermions, EE can be obtained from the cumulants of the number of transmitted 
charges distribution. However, this result may not hold in the general interacting 
case. 

3.2. Different kinds of entropy 

Without trying to be exhaustive, let us discuss a few different measures of entan- 
glement that are of interest for quantum spin systems, such as described by Eq. ([T]). 

3.2.1. Von Neumann entanglement entropy 

The von Neumann EE {S^^) quantifies the bipartite entanglement between two 
parts of a quantum systenJSl. For a quantum state \^), EE between a part and 
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the rest of a system is Sn = — Trpj^lnpji where pn = Trf2|\I')(5'| is the reduced 
density matrix of Q obtained by tracing out the rest of the system Q. An important 
property of EE is that Sq = Sq and it appears naturally that Sq is only related 
to the common property of fl and fi, their boundary: general arguments indeed 
indicate that Sq typically scales with the size of the boundary (so called area 
law)^Foi critical systems however, logarithmic corrections can be present. If the 
critical system is conformal invariant, the amplitude of the logarithmic corrections 
is related to the central charge of the corresponding conformal field theory (CFT). 
This has been shown in one dimension (IdjEEand for some CFT in two dimensions 
(2d;P, in which case the coefficient also depends on the geometry of O. The scaling 
of the EE with the size of O therefore contains precious informations on the state 
of the physical system and can be used e.g. to detect criticality. This has been 
successfully demonstrated for instance in quantum spin systemsEEEEH 

3.2.2. Renyi entropy 

In quantum information theory, the Renyi entropy is a generalization of S^^ for 
any integer n > 1: 

^n-^lnTrpJ^. (4) 

In fact, S^^ can be viewed as a limit of Sn when n — > 1. 

It turns out that Renyi entropy behaves quite similarly as S"^ . In particular, 
for critical Id systems, it has also logarithmic correction: for periodic boundary 
conditions (PBC) and a block of x sites, one haSH 

Sn = ^{l + -)lnx + K (5) 
6 n 

(with K a constant), thus allowing to extract the central charge c. 

Recently, a crucial point has been realized: it is possible to compute Sn for 
n > 1 with QMC in the VB basiJ^^. Despite a sampling problem that challenges 
the obtention of Renyi entropies for all blocks in the same simulation, the authors 
of Ref. [in] have provided compelling evidence in favor of a strict area law for 5*2 in 
the 2d Heisenberg case. 

3.2.3. Valence Bond entanglement entropy 

For a given bipartition of sites and bipartite VB state |<I>), let us consider a subsys- 
tem n. One can define the Valence Bond EE (5^^) of this state as: 

SX^(|$))=ln(2).na(|<i>)) 

where n^(|<I>)) is the number of singlets that cross the boundary of OlHEll The 
constant ln(2) is used to match the EE for a single site. In general, the GS of the 
Hamiltonian H is not a single VB state. For a linear combination I^P) = 
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with 1$,) bipartite VB states, we have defined 5^^(|*>) Y.^aiSn^{\^^))/Y.ia^, 
which can thus be computed from the sampHng of the GS of H by using QMCEH 
It turns out that S^^ has a similar behaviour as S"^^ for gapless and gapped 
Id systems (with respectively existence and absence of log corrections) . In 2d, S^^ 
satisfies a pure area law for gapped ground-state but has multiplicative log correc- 
tions in the gapless Neel phas^HHI that seems to differ from S^^ which has been 
claimed to satisfy a strict area law, according to DMRG simulations using ladders 
with many legs.^^ 

Recently, Lin and Sandvik have also proposed to measure a different VB EE by 
computing the average of the number of crossing links in the GS: 5'VB,2 ^ (rt^). 
This quantity scales as S'^^, but with smaller subleading scaling corrections for a 
spin chainj^fl. 



3.2.4. Loop EE 

For a single VB configuration, its loop graph is identical to the configuration itself. 
Therefore, EE is identically equal either to the number of bonds crossing the bound- 
ary, or to the number of loops crossing it. Using this idea, Lin and Sandvik have 
proposed to extend the definition to an arbitrary state by computing the average 
number of loops that cross a boundary between two blocks.!^ 

It turns out that in Id, the scaling is logarithmic with the size of the block with 
a prefactor ln(2)/3 for uniform chains with PBC, while in 2d, 5''°°? satisfies a strict 
area law for the Neel state.l^ 



3.2.5. Other proposals 

If one considers a conserved quantity (like the total number of particles in the 
canonical ensemble, or the total magnetization for a spin system), then it could be 
useful to investigate the fiuctuations of this quantity on a finite part of the systcm.t2il 
However, contrary to EE, one needs to make a particular choice that may depend 
on the physics so that this approach looks less systematic. 

Recently, it has also been proposed to look not only at the entropy, but at the 
whole distribution of eigenvalues of the reduced density matrix of a block, namely 
the entanglement spectrurrPH. Indeed, the entanglement spectrum of a block $7 
appears to be related to the physical spectrum of the same Hamiltonian if the system 
was restricted to this block. Because of that, it can be used to investigate edge 
physics of quantum Hall effect even if the total system has no edgeJ^H, Concerning 
magnetic systems, this strategy has been pursued on spin-1/2 ladders by measuring 
the entanglement spectrum of one leg, that shows Id-like features^'*. Unfortunately, 
because the full spectrum is needed, only ED techniques can be used, which does 
not allow to reach large systems. 
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3.3. Is it well defined ? 

While S^^ or Renyi entropies are clearly good quantum observables, and loop 
entropies are defined in terms of geometric objects (number of bonds or loops that 
cross a given boundary), but because of the overcompleteness of the VB basis, it 
is not granted that one would get the same results if the wavefunction was written 
down in a different (but equivalent) way. 



3.3.1. VB entropy 

At first glance, S^^ seems ill-defined as the bipartite VB basis is overcomplete. 
1^) can indeed be rewritten as as a different linear combination of other states 
for which could be different. However, turns out to be conserved in any 
linear combination between the bipartite VB states. It can be indeed showrP^ that 
for any linear relation J2i = between the bipartite VB states we have 

J2i (^i^n^ = for all fl. Another issue is the choice of the bipartite basis: in 
the general case, does depend on the precise choice of bipartition. However, if 
l^f) satisfies a Marshall sign criteriorPsl, a genuine bipartition where all Ui > (in 



the notations of Sec. 3.2.31 exists and should be taken. Not all singlet states satisfy 
this, but this is the case for e.g. the GS of Eq. ^ on bipartite lattices.'^ 

An alternative definition of the occupation number, which moreover is also valid 
for any spin S, has been given in Ref. 1271 Let us first define a spin-S" dimer as the 
two-sites singlet state: 

= E (-l)^--|-s„+s.) (7) 



We also define the reference state \Rs) = \ — S, +S, . . . , —S, +S) where the state 
is written in the (SiSz eigenstates basis and ordered such as A and B sites appear 
in alternating order. In particular |i?i/2) is nothing but the Neel state. As already 
noticecPfor S — 1/2, the reference state has an equal overlap with all bipartite VB 
states: (Rsl'Pv) does not depend on the bipartite dimer covering T). This property 
is established by a direct evaluation leading to (Rsly^v) = ^/{"^S + 1)^/"*, with N 
the total number of spins in the system. 

After some algebrgP^, one obtains that the "average" number of VBs between 
site i and j is given by: 

The total number of VBs crossing the boundary of a given bipartition Q, and there- 
fore iS*^^, are then easily obtained by summing over all i in f2 and j in f2. Note 
that this formulation is explicitly independent of the linear combination chosen to 
expand Ivf) on the overcomplete bipartite VB manifold as it only involves projec- 
tions of 1^). Moreover, this definition of S^^ as a projection is VB basis-free and 
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therefore aUows numerical computations outside the VB QMC scheme such as with 
ED and DMRG (see Sec.|4]). 

In summary, is a perfectly well-defined quantity for any singlet wave fonc- 
tion, although it does depend on the choice of bipartition. 

3.3.2. Loop entropy 

The previous discussion emphasizes the fact that any quantity designed for measur- 
ing entanglement should be checked first to be consistently defined with respect to 
quantum mechanics. In particular, for a given singlet state, it should not be tied to 
a particular VB representation among all the possible choices made possible by the 
massive overcompletness of the (bipartite) VB basis. As obvious as such a require- 
ment may seem, it is worth mentioning that since the various proposals have been 
introduced in the context of QMCliZli5122|, they are constructed as statistical aver- 
ages which do not necessarily make them properly defined as quantum mechanical 
averages (expectation value of an hermitian operator) or projections (like S^^). 

A simple argument based on a 6-sites example shows that S'loop does not sat- 
isfy the criterion raised in the previous paragraph. Despite its manifest physical 
contenlpSl, it is deeply tied to its stochastic definition inside the QMC scheme and 
cannot be consistently defined as a quantum average nor a projection. We con- 
sider two blocks A = (1,2,3,4) and B = (5,6) (see Fig. [T](a)) and the six bipar- 
tite states la), |/3), I7), \5), \e) and \C,) (see Fig. [TJb)). Due to overcompleteness, 
|(^) = |a) -f 1/3) + I7) and \^) — \5) + \e) + |C) are two representations of the same 
singlet state, \tf) = |-0) . Note that the example is intentionally built such as the 
weights involded in the decompositions \lp) and ji/;) are positive in order to stay 
closer to the situation encompassed by the projection QMC scheme. As shown in 
Fig. [TJc), using the notations of Ref. HOI {^\^abW) = 11/2 while (iy9|A^s|'0) = 6, 
excluding the possibility for 5'ioop = {^ab) to be a quantum mechanical average. 

Let us remark that even if this point raises a difficulty in the interpretation of 
•S'loop from a quantum mechanical point of view, this lessens in no way the numerical 
evidence that S'loop and S^^ probably have a very similar behavior. 

4. Numerical results 

As an illustration, we now present some numerical results for S"'^ and for 
frustrated spin-1/2 and spin-1 chains. 

4.1. Heisenberg spin-1/2 chain 

L 

Let us consider the frustrated spin-1/2 chain H = JiS^.S^+i + J2Si.Si-|_2 where 

2 = 1 

we set Ji = 1 and vary J2. The physics of this spin chain is well understood: for J2 
smaller than the critical value J| ~ 0.241167 (Ref. [55|), the system displays anti- 
ferromagnetic quasi-long range order, with algebraically decaying spin correlations. 
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Fig. 1. (Color online) 6-sites example showing that Sioop introduced in Ref. I20l is not independent 
of the chosen decomposion of a given singlet state on the overcomplete bipartite VB basis, (a) 
Definition of the subsystems A = (1,2,3,4) and B = (5,6). (b) The six bipartite VB states, (c) 
Loop contributions {v\Aab\^') = (number of crossing loops) X {v\v') (red numbers) involved in the 
computation of Sioop for a single singlet state using either the representation \ip) = |q) + \(5) + I7} 
or = \5) + \e) + Note that, due to overcompleteness, \ip) = \ip). 



For J2 > J2, the system is located in a gapped dimerized phase which sponta- 
neously breaks translation symmetry. In Ref. |271 we have studied both and 
S^^ in both phases using large-scale DMRG technique. 

To summarize our findings, we observe a similar behaviour for both EE in the 
gapless or gapped phases with respectively logarithmic increase or saturation. If 
we fit the prefactor of the logarithmic law in the gapless phase to a form S^^ = 
Ccff/31n(a:') -|- Ki and S^^ — ■jcS ln(a;') -t- K2, we obtain results very close to GET 
prediction^ZEEH c = 1 and 7 = 41n(2)/7r^. But even more interestingly is the fact 
that the values closest to these theoretical predictions are found to be precisely at 
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J|, which indicates that both quantities can detect the quantum phase transition. 

4.2. Bilinear-biquadratic spin 1 chain 

How can we extend the definition of the VB occupation number (and VB entropy) 
to spin values S > 1/2 1 Naively, one would try to count the number of 2-sites 
singlets formed by 2 spins S. While this is possible using the derivation outlined in 



Sec. 3.3.1 it is crucial to remark that, in contrast to the spin 1/2 case, such VBs 
do not form a basis of the singlet space of N spins S: in fact they form a basis 
of SU(N) singletd^, a symmetry which is only met by specific spin models with a 
symmetry higher than the usual SU(2). 

Here we propose an alternative, more generic, definition which consists in com- 
puting the "VB entropy of spins 1/2" in a spin S system. Without loss of generality, 
we will describe this method for systems with spin S* — 1. At each lattice site i, 
decompose the spin 1 S.^ into 2 spins 1/2 (denoted by i and 2), and consider 
the ensemble of VBs that can be formed by all these spins 1/2. At each site, the 
Hilbert space is of size 4 (instead of 3 for a single spin 1) and we will force the 
use of the physically relevant states by refusing that two spins 1 /2 form a singlet 
on the same lattice site. For systems on a bipartite lattice, this is easily done by 
taking the bipartite VB basis for the spins 1/2, where all A (B) spin 1 sites only 
host A (B) sites for the spins 1/2. All spin 1/2 VB states span the singlet sector 
of this fictitous spin 1/2 system, and therefore the full singlet sector of the spin 
1 system. We can now count the occupation of these spin-1/2 VBs using the for- 
mula given above. It is easy to see that the "spin 1/2 VB occupation number" 

S— 1 /2 

rtj j of a spin for two spins 1 at sites i (in sublattice A) and j (in sublattice B) 
in a wavefunction \'^), is obtained by putting together all spin-1/2 s operators : 
4,iS~i + 4,iS~2 + 4,2sJa + «t2Sj~2' resulting in 



s-1/2 (■Ril'S'+^rl* 



{Rim ■ 



Here \Ri) is the reference state for S = 1 defined in Sec. 3.3.1 Note the resemblance 
with the counting of spin-1 VB, except for the factor 2. Indeed, there is at most 
one spin-1 VB between two spin 1 sites, whereas there can be up to two spin-1/2 
VBs. With this definition, we can now compute the "VB entropy of spins 1/2" 
by counting how many spins 1/2 VBs emerge from a given block size of a spin 1 
system, multiplied by ln(2). It is important to remark that the spin 1/2 artificial 
system is useful only to identify our computations, but that all the simulations and 
numerical evaluations can be made directly in the spin-1 model and formalism. 
For illustration, we consider the spin-1 chain with bilinear and biquadratic terms 

•H = ^ COS0 S,.S,+i + sine (S,.S,+i)", (10) 

i 

which phase diagram is well known thanks to analytical and numerical studies!^ 
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Here, we focus on one side of the gapped Haldane phase, for instance by choosing 
6* = (Heisenberg point), — 10° and 9 — arctan(l/3) which corresponds to the 
famous AKLT point with an exact ground-statd^ and a short correlation length 
^ = l/ln3. On the other side, we want to see whether S^^ can detect the phase 
transition that occurs at 6 ~ — vr/4, the Babujian-Takhtajan (BT) point, where 
criticahty is known to be described by a SU{2)2 CFT with central charge c = 3/2!^ 
In what follows we use DMRG to compute the ground-state of Ti. and the von 
Neumann and spin-1/2 VB entropies for chains of size L = 64 with PBC. 




1 ' 1 ' 1 ' 


' 1 ' 1 ' 1 


— o 
o 

~o 


—-^^—^•"^ - 

o — 
o 

6~ 


^ 


^ 




or 




e=io° - 


-o 

1 1 1 1 1 1 


□- 

1 1 1 1 1 1 



10 20 30 40 50 60 



2 

1.8 
1.6 
1.4 
1.2 
1 



2.5 



1.5 



1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 


□ 
□ 

□ 

□ 


□ 
□ 

□ 

□ 


_0 




o_ 


6> 




Heisenberg 




P , 1 , 1 , 1 , 1 


,1,1° 



10 20 30 40 50 60 




d(x)=L/7tsin(xJl/L) 



Fig. 2. Scaling as a function of block size x of the von Neumann (black circles) and spin- 1/2 VB (red 
squares) entropies for the bilinear-biquadratic spin 1 chain for various values ol9 : d = arctan(l/3) 
(AKLT point), B = 10°, 6 = (Heisenberg point), and 6 = -tt/4 (BT point). In the last case, 
entropies are plotted versus conformal distance d{x) = L / tt sm(x-n / L) on a semi- log scale and lines 
are best fits to a logarithmic divergence (see text). We used L = 64 and PBC. 



At the AKLT point, S*^^ saturates to 2 In 2 rapidly (with the same correlation 
length ^ as the bulk one) as was known from previous study.El The VB entropy of 
spins 1/2 is constant for all block sizes at the AKLT point, and the constant is also 
exactly equal to 21n(2): this is easily understood noting the VB solid structure in 
the AKLT ground-state where there is exactly one spin 1/2 VB on every lattice 
link. For the Heisenberg point or 6 — 10°, saturations of both entropies can be 
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readily seen too, but now there is no direct relation between the constants towards 
wich both S and S^^ saturate. 

For the BT point, S^^ is expected to show logarithmic behaviour with a slope 
c/3 (since we are using PBC), with c = 3/2 from the known criticality. Indeed, we 
observe a very good agreement with our data. The VB entropy of spins 1/2 also 
exhibits a logarithmic divergence with block size, with a slighlty different prefactor 
75^(2)2, exactly as in the spin-1/2 SU{2) case. Theoretical work by Saleur and 
JacobserPSl, made along the lines of Ref. predict Jsu{2)2 ~ 2. 7 = 81n(2)/7r^ ^ 

0. 562. again close but different from c/3 — 0.5. Our best fit gives 75(7(2)2 = 0.667, 
slightly different from the theoretical value. This discrepancy is likely due to log- 
arithmic corrections due to the enhanced symmetry at this point, as was found in 
the spin 1/2 Heisenberg casd^. 

5. Conclusion 

It is now obvious that entanglement measurements can shed light on the nature 
of the physical state, without having any a priori input. While the computation of 
^vN \g^j.gQ interacting 2d systems remain a very difficult task, there has been 
several recent progress using large-scale QMC simulations of spin systems in the 
VB basis. We have reviewed some proposals and discussed especially if these new 
estimators are valid quantum observables or projections. We focused on the VB 
EE, which can be either defined as a geometrical quantity, or as a projection onto 
a particular state, and which has been shown to exhibit several interesting features 
for gapped or gapless spin-1/2 phases. We have described how it can be extended 
to arbitrary spin and provide numerical data for frustrated spin-1 chain, where a 
critical state can be detected using this tool. Clearly, more work should be devoted 
to understand similarities and differences between these various estimators. 
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